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Abstract 


The effect of thin turbulent layers, including boundary layers 
and shear layers, on light propagation is examined from a theoretical 
point of view. In particular, a mathematical model is developed to 
describe the interaction between the aerodynamic or, more precisely, 
the density fluctuations and the electromagnetic field. It is 
assumed that the turbulence induces a normally distributed phase 
aberration which is a homogeneous random function in the plane of the 
aperture. Hypotheses concerning the density fluctuations in the 
layer sufficient to guarantee such a phase aberration are exhibited. 

The optical degradation is described in terms of the optical 
transfer function (OTF) and the Strehl ratio (I/I Q ) which are random. 
Expressions for the first and second moments of these two parameters 
are developed from the definitions. Asymptotic ("large" aperture) 
approximations to these expressions are derived and discussed. 
Finally, the exact and approximate results are compared for several 
"typical" values of the ratios of aperture diameter to scale of 
density fluctuations and rms phase aberration to wave length 
respectively. 
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1. Introduction 

The effect of turbulent layers on the propagation of electro- 
magnetic waves has drawn much attention over the past three decades. 
In 1950, Booker and Gordon [1] developed a theory which explained 
the scattering of radio waves in the troposphere. In 1956, Stine 
and Winovich [2] conducted an experimental investigation of light 
diffusion through turbulent boundary layers which tended to validate 
the applicability of the Booker-Gordon theory in this new context. 
Then, in 1971, Sutton [3] analyzed the Stine and Winovich data 
further and concluded that, for imaging apertures much larger than 
the turbulence scale size, scattering decreases contrast but not 
resolution. 

While the above investigators were examining the influence of 
turbulence on resolution by studying the induced scattering, others 
looked at the influence on the optical (or modulation) transfer 
function. Hufnagel and Stanley [4] developed an expression for the 
average or expected value of the MTF . In particular, they showed 
that the average transfer function can be written as the product of 
the diffraction limited transfer function and an attenuating factor 
which incorporates properties of the turbulence. Fried elaborated 
on this theme in a series of papers [5, 6,7,8] in the late 1960's, 

Much of the recent theoretical development connecting optical and 
turbulent parameters has been summarized by Wolters [9] . 

In the present work, we will attempt to present the fundamental 
theory relating turbulence statistics and mean values of the optical 
parameters, describing in some detail the relevant statistical 
assumptions. Also, we will discuss the second moments of the optical 
parameters when the aperture is much larger than the scale of 
turbulence . 
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2, Mathematical Model of Aero-Optic Interaction 


The model described below does not originate with the author, 
but appears, at least implicitly, in much of the current literature 
on aero-optic interaction. Perhaps the most complete discussion has 
been given by Wolters [9] . However, because the assumptions and 
hypotheses aren't always explicit, and because the development is 
not easily accessible in the literature, it is believed that the 

somewhat more thorough discussion here is warranted. 

We begin by assuming that the turbulent layer of thickness L 
induces a phase aberration but no amplitude degradation of the optical 
wave front (monochromatic of wave length A) . Hence the wave field 
u on the aperture of the receiving optics* can be described as 
follows (see Figure 1) : 

u(x,y) = jpxp [ikA (x,y)], (2.1) 

A(x,y) =jn(x,y,3)ds 

= ju+n 1 (x,y ,a) ] da 


where 


= L + 


n 1 (x,y,a) da 


= L + A 1 (x,y) . 


In the above, k is the wave number (2 w/a) , n is the local refractive 
index, and n^ is the deviation of the refractive index from its value 


* 

Even though we shall be discussing a passive receiving system, 
it should be noted that an active or propagating system (e.g., 
a high energy laser beam) would experience similar degradation. 


256 




in a vacuum. The expression (2.1) is the first order geometrical 
optics approximation, where it has been assumed that |n^|<< 1 [10]. 

The phase shift A is sometimes called the optical length of the 
ray's path through the turbulent layer. (See, for example, [ll,pll5].) 

Without making some assumptions regarding the statistical nature 
of the random phase A^, further progress would be difficult if not 
impossible. In particular, we make two assumptions: 

(i) The random variable A^(x,y) is normally distributed 
for any x and y, 

and 

(ii) A 1 (x,y) is a weakly homogeneous random process; 
that is, 

<A^(x,y) > = constant (2.2) 

<(A 1 (x 1 ,y 1 ) - <A 1 >) (A 1 (x 2 ,y 2 ) - <& 1 >)> = R(x 2 ~x ;l ,y 2 ~ y-j^) , 

where < > denotes ensemble average or mathematical 
expectation. 

Assumption (i) can be justified by considering the relationship 
between the optical parameter A^ and the air density p . By the 
Glads tone-Dale law we have 

n-L^y,*) = G p (x,y,«) , (2.3) 

3 

where G = 0.000223 m /kg is the so-called Glads tone-Dale constant. 
Hence, substitution of (2. 3) into (2.1) yields: 
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A 1 (x / y) = G 


(2.4) 


L 

J p (x,y,a) da. 

0 

Assume that L >> £ a , where £ g is the integral scale of the density 
fluctuations in the direction of wave propagation. It follows 
that we can partition the random layer into strata of thickness h, 
with h (and hence N) chosen so that 


£ < h 

a 

We then have that 

A x (x,y) = G 


= L/N 

<< 

L . 

N 

C 1 


jh 

> 


p (x,y,a) da 

Z 

— J 



1 >- j 

j=l (j-l)h 


(2.5) 


But the random variables determined by the integrals in the sum 
(2.5) are approximately mutually independent (since h > £ a ) ; then, 
since N >> 1, the Central Limit Theorem [12, p.266] gives the 
desired result. It should be observed that the assumption L >> £^ 
is often satisfied; indeed, in practice, 

10 < L/£ < 40. 

2 

(See, for example, [13 , Tables 1 and 2].) 

The first of conditions (2.2) follows from (2.4) if we assume 
that < p(x,y,a) > depends solely on a and not on the aperture 
coordinates (x,y) . The second condition, translation invariance of 
R, follows from the assumption that the covariance of the density 
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( 2 . 6 ) 


fluctuations, R , satisfies 


R P (X 1^1^1 ;X 2^2' H 2 ) = R p (x 2" x l'y2" y l ; * 1 '* 2 ) 


Then 


r(x 2 -x 1 , y 2 ~yi) = G ' 


L L 
r 


R (x 2 -x 1 , y 2 -y 1 ; a 1 ,« 2 )d« d« 2 . (2.7) 


0 0 


The assumption (2.6) is consistent with experimental evidence which 
suggests that the rms density and the turbulence scale length in 
any direction (i.e., * x , % , or £ g ) vary with a [13] . 

In practice, it has been found that R p can, with good accuracy, 
be represented as follows [14, Figure 20J : 


R p (x 2~ X l' y 2 -y l ; S 1 ,S 2 ) = R p (x 2 _X l ,y 2~ y l ? U/V) 


( 2 . 8 ) 


= ct (u) exp \ - 


x 2 -x l 1 +/ y 2 y l 

a y (u) ' 


w 


+ 


a z (u)i 


where 


u = (■z 1 +a 2 )/2 


(2.9) 


v = •z 2 - Bl . 


But, making the change of variables (2.9) in (2.7) we find 
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( L/2 2u 

r r 


R(x 2 -x 1/y2 - yi ) = 


G 2 < 


R p (x 2" x l' y 2“ y l'* U/V) dv 


du 


0 -2u 


V L 


2 (L-u) 


R p (x 2 -x 1/y2 - yi ; u,v) dv du 


> . ( 2 . 10 ) 


L/2 -2 (L-u) 


The relations (2.8) and (2.10) allow calculation of R once the -z 
2 

variations of a .£ , £ , and £ have been determined. Fmall y , 

p x y s 

from (2.8) and (2.10) it can be shown that 


a = R(0,0) * 2G' 


a (z) £ (z) dz , 

P Z 


( 2 . 11 ) 


for L > > £ 

2 
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3. 


Second Moments of Optical Parameters 


We proceed now to develop expressions for the second order 
statistics of the optical transfer function t. We begin with the 
definition : 


r 


x (x,y) = ( 4/ir) 


u(5+ n + u*(£- n- ^-) d£dn. 


(3.1) 


x f y 


where A is the area common to two identical apertures (of unit 

x / y 

diameter) displaced relative to each other a distance x and y along 
the 5 and n axes respectively. Note that t (x,y) = 0 if x^+y^ _> 1. 
Clearly, since u is random, then x is itself random. Hence, taking 
the expectation of (3.1) we obtain: 

<x (x,y) > = (4/rr) 




<u(£ + j, n+ ^) u*(£- n- ¥?) > dg dn 
x,y 


= (4/ it) 


= (4/ir) 


< exp ik [A 1 (?+ j, n+ j) -A 1 ( n - ^-) ] j > d? dn 


X/Y 


exp -k 2 [R ( 0 , 0) - R(x,y) ] d? dn 


x,y 


= exp j -k 2 a 2 [1 - r (x,y ) ] j x o (x,y). 


(3.2) 


2 

where ct = R(0,0) is the variance and r(x,y) is the correlation 
coefficient of the phase aberration A^ and t q ( x,y) is the diffraction 


* 

The aperture coordinates 5,n ,x, and y have been normalized with 
respect to the aperture diameter, D. Then (x,y) correspond to spatial 
frequencies (f x ,f y ) = (Dx/AR, Dy/AR) , where R is the focal length 
of the optical system. 
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limited optical transfer function associated with the receiving 
optics. We have invoked both assumptions (i) and (ii) of Section 2 
in obtaining (3.2). Hence, the average optical transfer function is 
just the diffraction limited optical transfer function t q attenuated 
by the factor 


T att (x,y) = exp } - (2ira/A) “ [1 - r(x,y)] ] . (3.3) 

In general, to describe completely the random process x(x,y), 
we must obtain correlations of all orders in addition to the 
expectation given by (3.2). This has in fact been accomplished by 
Barakat [ 15] . We will be concerned here however only with the second 
order correlation or auto-covariance function. But by following 
the procedure outlined in the development of (3.2), it is a 
straightforward task to verify that the auto-covariance function for 
x is given by: 


R 

T 


(x,y,x' ,y ' ) 


t (x,y) t (x' ,y ' )(4 /tt) 
att att 


2 


rr rr 


F(u, v) 


yy yy 




d£dnd£ ' dn ' , 


(3.4) 


where 

u = C'- K 
v = n n 


and 
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I 

F(u,v) = exp I -(2t7 ct/X) 

-r(u+ v+ ^ ) -r (u- v- Z-JX ) 

«<- 2^2, V- ^ )] 

It should be observed that since the covariance function R^ depends 
explicitly on x,y,x' , and y' rather than on the differences x'-x and 
y'-y, the optical transfer function t is not a homogenous random 
process. Indeed, since we are considering only phase aberrations, 
lx(x,y) | _< T Q (x,y) _< 1 . Hence t is clearly not even normally 
distributed. A fuller discussion of this matter can be found in (15], 
Before proceeding, note that the variance of t ( x,y) is given by 

a^(x,y) = R T (x,y;x,y) . (3.5) 

We turn now to the so-called Strehl ratio, I/I Q , the ratio of 
the maximum intensity in the image plane with aberrations, I, to the 
maximum intensity without aberrations, I Q , both in response to a 
point source- We will use the following definition (See, for example, 
[ 16 , P 88] ) : 

oo oo 



r (u+ 


x +x 




x (x,y) dx dy / 


t q (x,y) dx dy 


— oo —c>o 

Now since T Q (x,y) is defined as the normallized convolution of the 
unaberrated pupil function, it can easily be shown (see, for example. 


[17 , p 166] ) that 
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oo 
r r 


T 0 (x,y) dx dy = tr/4 . 


- oo 

Hence (3,6) can be rewritten as: 

oo 


I/I o = (4/ it) 


t (x,y) dx dy . 


(3.7) 


— oo 


It follows easily from (3.7) that the first and second moments of 


the Strehl ratio are given by: 

» 

<x(x,y) > dx dy 


< I/I Q > = (4/tt) 


OO 
” r 


(3.8) 


— OO 


and 


oo 

err. 


a2 I/I = (4 / 1r) 2 


— OO 


(x,y,x' ,y ) dx dy dx' dy' 


(3.9) 


respectively. It is interesting to observe that if the fluctuations 

in t at different spatial frequencies (x,y) are perfectly correlated, 

then R t ( x,y,x' ,y ' ) = a T (x,y) c t ( x' ,y ’ ) and hence, from (3.9), 

OO 


a 


I/I 
' o 


(4/tt) 


ct (x,y) dx dy 

T 


— oo 


(3.10) 


In this case, the rms Strehl ratio is determined simply by the volume 
under the rms optical transfer function. However, in general, the 
optical transfer function does not fluctuate uniformly over its width 
and so (3.10) cannot be expected to give accurate results. 
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The equations for the moments of t and I/I developed above 
could in principle be used to obtain quantitative results for a 
given turbulent layer by utilizing a high speed digital computer to 
ca,rry out the required integrations. This procedure can be quite time 
consuming however, especially for the four dimensional integrations 
of (3.4) and (3.9). However, by pursuing an asymptotic (D f =>) 
analysis, an approximation to these integrals can be achieved. This 
approach will be pursued in Section 4 below. 

Before proceeding, we should remark that it is of some interest 
to consider the modulus of the optical transfer function, |x|, sometimes 
referred to as the modulation transfer function or MTF. Since x is 
in general a complex- valued function, the MTF does not yield information 
concerning the phase of x . However, in recent exneriments , Kelsall 
has utilized the fast shearing interferometer [18] which measures the 
MTF rather than the desired optical transfer function. Although the 
two parameters are obviously not equivalent, there is a relationship 
between the two. In particular, it is not difficult to show from 
their respective definitions that: 

< x > < < MTF > < x (3.11) 

— — o 

and 

2 2 2 , .. 2. 2 

CT < cr < x (1 — X. )< X . (3.12) 

MTF T ° att ° 

Clearly, from (3.11) and (3.12), the MTF suffers less degradation 
than does the optical transfer function x . Note that even though, 
from (3.2) , < x > is real, this does not necessarily imply that 

< x > = < MTF > . 
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4. Asymptotic (Large Aperture) Approximations 


The equations developed in Section 3 allow one, in principle, 
to calculate the first and second moments of the optical transfer 
function, t, and the Strehl ratio, 1/I C . Indeed the expression (3.2) 
for <x> is explicit and needs no further comment. However, the 

2 2 

calculation of a T , R t , <1/1,-, > r and require multiple numerical 

integrations which can consume large amounts of computer time. 
Furthermore, the exact expressions tend to conceal the influence 
of variables like rms aberration, a, and aperture size, D, on the 
parameter in question. Hence, in this section, we will develop 
approximations for the moment expressions derived in Section 3 for 
the case D>>s, . 

A 

Consider first the average Strehl ratio? by (3.2) and (3.8), 
we have : 


oo 


< I/I Q > = (4/ir) \ \ exp ( -k a [l-r(x,y) ] 


.2 2 


T Q (x,y) dxdy 


- OO 

= (4/tt) exp (-k 2 a 2 ) 


OO 


exp j k 2 a 2 r(x,y) 


T o (x,y)dxdy (4.1) 


— oo 


Now expand the exponential part of the integrand in a Taylor series, 

CO 


2 2 

exp Ik a r(x,y ) 


I 

p=0 


n 2 2 \ 
(k a ) 

~pT 


r F (x,y) . 


Then, substitution into (4.1) yields: 
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< I/I c > 


= exp (-k 2 cr 2 ) ( 


oo 


1 + 


P=1 


(k 2 a 2 )P 


OO 
r 


r? {x,y) T Q (x,y) dx dy 




— oo 


(4.2) 


The expression (4.2) is exact; no approximations have yet been made. 
Now let's assume that r is of the form* 


r (x,y ) = exp 


- ~\[( D/* x ) 2 x 2 + (D/l ) 2 y 2 


(4.3) 


We will consider the limit D/& x + 08 while & x /£y remains constant. 
Then (4.3) can be rewritten 


r(x,y) = exp 


- (D / V “\/ x2 + U x /A y ) 2 y 2 


(4.4) 


Clearly, as D/fc grows larger, the graph of r(x,y) becomes narrower. 
In fact, it is an easy matter to verify (keeping in mind that the 
volume under the delta function 6 (x,y) is unity) that 


rP(x,y) v 


2t T 


P (d/x, x ) (D/a y ) 


6 (x,y) 


(4.5) 


as I 00 / for any p > 0. The details of derivation of (4.5) can be 

found in the appendix. 


Substitution of (4.5) into (4.2) gives 


2 2 , 


< I/I > 'v, exp (-k a ) / 


1 + 


oo 

S’ ' ,, 2 2. p 

x (k a )^ 


oo 

r 


2tt 


, Pl P 2 (D/A x ) (D/£ v ) J 
p=i x y 


6 (x,y) t Q (x,y) 


— oo 


This requires either that the scale lengths and & be constant or 
that ^ x (u) and £ y (u) be replaced by average values in (2.8). 
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= exp (-k 2 a 2 ) ( 1 + 


8 


(D/JU(D/*) p 2 pl 

x y p=i 


oo 


I 


(k 2 a 2 ) P 


(4.6) 


as D/^ x t <*>. Since the sum occurring in (4.6) will appear again. 


let us define 


oo 


F(x) = 


P=1 


2 , 
p pi 


(4.7) 


(Note that this sum is convergent for all x.) Then, from (4.6), we 
have / \ 


< I/I Q > ^ exp (-k 2 a 2 ) { 


i I 8 , , / 1 2 2 . 

1 WCTWTT F(k ° > 

y 




(4.8) 


as D/ ' Z x + ”* ( For reference, the function F has been graphed over an 

interval of x sufficient for most conceivable aberrations; see 

Figure 2.) It might be observed that the second term of (4.8) can be 

viewed as the contribution from what I-Iogge and his colleagues called 

the incoherent beam [19] . They concluded that for a phase-aberrated 

beam, the far-field irradiance distribution can be written as the 

sum of two beams; one beam is the diffraction limited beam attenuated 

2 2 

Joy the factor exp (-k a ) and the other beam is much wider and 
contributes an amount, on-axis, proportional to the second term of 
(4.8) . 

We turn our attention now to the variance of the optical 


a ^(x,y) = t 2 (x,y) (4 /tt) 2 
att 


£ 

r 





exp -k 2 a 2 g(£,n,£' ,n') 

-l] 

J > 

J , 

) \ L 

1 


A A 

x,y x,y 


(4.9) 


268 



where 


gU, n, V , n') = r(u+x,v+y) - 2r(u,v) + r(u-x,v-y) 
u = 5'-5 

V = 1} '-'T) • 

Following the procedure adopted in the analysis of < 1/ I > t we expand 
the exponential in the integrand in a Taylor series and (4.9) becomes: 


^ T (x,y) = ^ tt ( x 


OO /■/* 

, y) g 


g p ( ?, n, S' ,n' ) d n dg' 


P =1 


A x,y A x,y 


(4.10) 


The expression (4.10) is exact but not very useful as it stands. It 
remains to estimate the 4-dimensional integral in the case 

O >> £ x • 

Again we assume that the correlation function r is of the form 
(4.3) and note that it approximates a delta function as D/Jl + <*> . 
Specifically, we have the asymptotic approximation given by (4.5). 
Hence, each of the three terms in g (defined by (4.9) approach a 
delta function in shape. But, if (x,y) f (0,0), each delta function 
is centered at a different point in the u-v plane. In this case, 
cross terms in the product g P can be neglected and we have 


g p u ,n , 5 ' ,n ') sr r p (u+x, v+y) + (-1) p 2 p r p (u, v) + r p (u-x, v-y) 


2tt 



P 2 (D/A x ) (D/* y ) 


$ (u+x,v+y) + (~1) P 2 P 6 (u,v) + S (u-x,v 


(4.11) 
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as D/2,^ + °° . Note that this approximation is least accurate near 
(x,y) = (0,0). Substitution of (4.11) into (4.10) gives 

(~l) p (k 2 o 2 ) p 


oo 


a 2 (x,y) x 2 (x,y) 

att 


2tt 

tT/ wywr 


r I 

P = 1 


2 , 
P Pi 


(4.12) 



<5 (u+x, v+y) + (-1) P 2 P S (u, v) +6 (u-x, v-y)j d§dnd£ 1 dn ' 


A x,y A x,y 
It remains to evaluate the integral. 

We proceed now to analyze the first integral of (4.12). The 
remaining two will follow directly. We have 



— oo 


5 (u+x, v+y) dSdnd^'dn' 


oo 

r 


j . 

~ OO 


G o (5+ f' n+ J )G o (5 “ 7 fT) ~ 7 > G o U_ ^jr) d?dn , (4.13) 


where G q is the unaberrated pupil function; i.e. 


G 0 (x,y) = 


1 / 

0, 


*^x 2 +y 2 < 1/2 


V 2" 2 

v x +y > 


+y > 1/2 . 
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From Figure 3, it can be seen that (4.13) is just the area of 
intersection of three circles of unit diameter. In fact, the middle 
circle contributes nothing, and we can rewrite (4.13) as follows: 


oo 


rr 


<S (u+x, v+y) d5drid£ ' dn ' = 


A x,y A x,y 


G U+ n + $)G ^,r\- ^) d^dn 


— oo 


= \ T o ( 2x, 2y) . 


(4.14) 


The remaining integrals in (4.12) can be evaluated similarly 
and we find: 


a (x,y) -V 


8 


T 2 (x,y) 


(D/£ x ) (D/2 y ) att 


F(2k 2 a 2 )T Q (x,y)+2F(-k 2 a 2 )T o (2x,2y) 


(4.15) 


as D/2. + <=° , where F is given by (4.7) . 

X 

One could conceivably carry out analogous arguments to estimate 
2 

and hence a . However this was not attempted. 
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5. Conclusions 


In order to judge the accuracy of the asymptotic approximations 
derived in Section 4, those expressions and the corresponding exact 
results from Section 3 have been calculated for several values of 
phase aberration a/X and aperture size D /^ x * These calculations were 
carried out on a Burroughs 6800 digital computer. 

Figure 4 illustrates the variation of average Strehl ratio with 
D/fl, for three phase aberrations. The exact Strehl ratio was 
calculated from (3.8) and (3.2) using, successively, Simpson's rule 
and the four-point Gaussian quadrature formula to evaluate the double 
integral. Note that if D << Z , the optics are essentially insensitive 
to the turbulent layer. Then, as D/l grows larger, the average 

X 

Strehl ratio decreases to an asymptote determined by the aberration 
a/ A. It is clear from Figure 4 and (4.8) that the so-called 
"infinite aperture" (or zeroth order) approximation 

< I/I o > ^ exp (-k 2<j2) 

is reasonably accurate for D > 6& x . The first order correction 
(given by the second term of (4.8)) varies from about 3% for a/X = 0.08 
to 22% for a/X = 0.2 when D = 6£ . 

It is appropriate here to relate the expression for < I/I 0 > 
given by (4.8) to the work of Hogge, Butts, and Burlakoff [19] . In 
obtaining (4.8) from (4.1), the exponential term exp [k a r(x,y)] 
was expanded in a power series. In [19], only the first two terms of 
this series were retained, thus limiting the validity of those results 
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to a/ A no larger than about 0.1. No such limitation applies to 
the results derived here. On the other hand, Hogge, et al, reached 
important conclusions regarding the shape of the focal plane 
irradiance distribution (discussed in Section 4) , whereas the 
spatial distribution of irradiance has not been considered here. 

Figures 5, 6, and 7 illustrate the correspondence between the 
exact (Equation (3.5)) and the approximate (Equation (4.15)) values 
of for successively larger D /*- x (3,6, and 10) with a/\ = 0.2. 

The exact expression (3.5) was evaluated using a Monte Carlo 
technique to approximate the four- dimensional integration. As 
expected, agreement is best for D/£ x =10. In fact, for this case, 
the accuracy of the Monte Carlo method employed to evaluate (3.5) 
is questionable and, hence, given the inherent computational errors, 
the exact and approximate values of can be said to agree. 

It is clear, especially from Figure 5, that 

c t (approx.) + “ as x/D + 0. 

This anomaly is the result of the assumption that (x,y) is not near 
(0,0) which leads to (4.11). 

The work described here suggests further research toward 
understanding the second-order statistics of the optical parameters 
I/I Q and x, namely: 

1. It is desirable to obtain approximations for and • 

Since the exact expression for these parameters involve 
complicated four-dimensional integrations , an approximate 
(closed form) expression would be especially useful. An 
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analysis similar to that described in this paper might 
prove fruitful. 

2. Perhaps by obtaining higher order asymptotic (D/& t «>) 

A 

approximations to < I /I > and better accuracy can be 
achieved for smaller values of D/Jl^. This would require a 
more accurate description of the aberration correlation 
r(x,y) than the 6-function analysis provides. 


3. 


Finally, the experimental data tends to substantiate the 
mathematical model employed here. For example, see 
Reference 13 for a discussion of the aperture scaling 
inferred by Figure 4. Also, a comparison can be made between 
the data and the expression for. < x > given by (3.2). 

However, there has been very limited effort expended to 


compute a . 


from measured data and virtually no attempt 


to compute c . Since these calculations can be accomplished 
routinely by modifying existing data reduction codes, it is 


strongly urged that this information be provided in future 
reports of experimental data. 
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APPENDIX 


The defining properties of the two-dimensional delta function, 

6 (x-x , y-y ), are as follows: 
o o 

oo 

c c 

f(x,y) 8 (x-x q , y-y o ) dx dy = f ( x 0 /Y 0 ) 

— OO 

8 (x-x q/ y-y Q ) = 0, (x,y) ^ (x 0/ y 0 ) . 

Although no ordinary function can satisfy these requirements, it is 

possible to construct sequences of functions, S^(x-x q , y-y Q ) , which 

approach the symbolic "function" 8 (x-x q , y-y Q ) as k + 

In particular, if 

oo 

c c 

s k (x_x o' y-y 0 ) dx d y = !» 

✓ J 
— CO 

for every k, and 


lim S k (x-x Q , y-y Q ) = 0, (x,y) ^ ( X 0 /Y Q ) / 

kt<=° 

then we say that 

lim S k (x-x o , y-y Q ) = 6 (x-x q , y-y Q ) . 

k + co 

Now, consider the function 

r P (x,y) = exp | - p"^(D/A x ) 2 x 2 + (D/£ y ) 2 y2 
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Then, after transforming from rectangular coordinates (x,y) to elliptic 
coordinates (r,0), where 


r 2 = (D/£ x ) 2 x 2 + (D/A y ) 2 y 2 , 


we have 


oo 

rr 


2ir 

r 


oo 

r 


r p (x,y) dx dy = 


- oo 


de 


(D/£ x ) 2 cos 2 e+(D/£ y ) 2 sin 2 e 


e" pr r dr 


2ir 


P (D/* x ) (D/A y ) 


Then, it follows that 


lixn P (D/ & x) <D/y 
2ir 

(d/a x ) +« 


r p (x,y) = 6 (x,y) 
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u (x,y) = 1 (PLANE WAVE) 



Figure 1. Phase-aberrated plane wave. 


4 Squares 
to the Inch 



280 








.8 


.6 


t ( x/D, 0) 


.4 


. 2 


0 



0 .2 .4 .6 .8 1.0 


x/D 


Figure 5. Second-order statistics of the optical transfer 
function for D/Jt =3 (Z =2,. a/\=.2). 




Figure 6. Second-order statistics of the optical transfer 
function for D/£ =6 (£. /£. = 2, a/X=.2). 


.8 


T (x/D , 0) 



x/D 


Figure 7. Second-order statistics of the optical transfer 
function for D/£ =10 (£ /£ =2, a/A=.2). 


